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Supernovae and Neutrinos 

John F. Beacom ^ 

'^NASA/Fermilab Astrophysics Center, Fermilab, Batavia, lUinois 60510-0500, USA 
beacomSf nal . gov 

A long-standing problem in supernova physics is how to measure the total energy and temperature of v^, Ut, 
V/j,, and i^T- While of the highest importance, this is very difficult because these flavors only have neutral-current 
detector interactions. We propose that neutrino-proton elastic scattering, u + p + p, can be used for the 
detection of supernova neutrinos in scintillator detectors. It should be emphasized immediately that the dominant 
signal is on free protons. Though the proton recoil kinetic energy spectrum is soft, with Tp ~ 2El/Mp, and the 
scintillation light output from slow, heavily ionizing protons is quenched, the yield above a realistic threshold is 
nearly as large as that from v^+p ^ + n. In addition, the measured proton spectrum is related to the incident 
neutrino spectrum. The ability to detect this signal would give detectors like KamLAND and Borexino a crucial 
and unique role in the quest to detect supernova neutrin os. These results are now published: J. F. Beacom, 
W. M. Farr and P. Vogel, Phys. Rev. D 66, 033001 (2002) | |arXiv:hep-ph /020522C| 1 ; the details are given there jl|. 



1. Introduction 

When the next Galactic supernova occurs, ap- 
proximately 10^ detected neutrino events are ex- 
pected among the several detectors around the 
world. It is widely believed that these 10^ events 
will provide important clues to the astrophysics 
of the supernova as well as the properties of the 
neutrinos themselves. Interestingly, recent break- 
throughs in understanding solar and atmospheric 
neutrinos each occurred when the accumulated 
samples of detected events first exceeded 10^. 

But will we have enough information to study 
the supernova neutrino signal in detail? Almost 
all of the detected events will be charged-current 
Ve + P ^ + n, which will be well-measured, 
both because of the large yield and because the 
measured positron spectrum is closely related to 
the neutrino spectrum. Because of the charged- 



lepton thresholds, the flavors 



and 



can only be detected in neutral-current reactions, 
of which the total yield is expected to be ap- 
proximately 10'^ events. However, in general one 
cannot measure the neutrino energy in neutral- 
current reactions. This talk presents an excep- 
tion. These four flavors are expected to carry 
away about 2 /3 of the supernova binding energy, 
and are expected to have a higher temperature 



than Vf, or Vf,. However, there is no experimental 
basis for these statements, and present numerical 
models of supernovae cannot definitively address 
these issues either. If there is no spectral signa- 
ture for the neutral-current detection reactions, 
then neither the total energy carried by these fla- 
vors nor their temperature can be separately de- 
termined from the detected number of events. 

But it is crucial that these quantities be mea- 
sured. Both are needed for comparison to numeri- 
cal supernova models. The total energy is needed 
to determine the mass of the neutron star, and 
the temperature is needed for studies of neutrino 
oscillations. At present, such studies would suffer 
from the need to make model-dependent assump- 
tions. This problem has long been known, but 
perhaps not widely enough appreciated. In this 
talk, I clarify this problem, and provide a realistic 
solution that can be implemented in two detec- 
tors, KamLAND (already operating) and Borex- 
ino (to be operating soon) . The solution is based 
on neutrino-proton elastic scattering, which has 
never before been shown to be a realistic detec- 
tion channel for low-energy neutrinos. 

In this talk, I will focus on just the problem 
of measuring the temperature and total energy 
of i/^, v-r, i^^, and P^., since everything else in 
understanding supernova neutrinos depends on it. 
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2. Cross Section 



3. Supernova Neutrinos 



The cross section for neutrino-proton elastic 
scattering is an important prediction Q of the 
Standard Model, and it has been confirmed by ex- 
tensive measurements at GeV energies (see, e.g., 
Ref. ||]). At the energies considered here, the full 
cross section formula [§||J^] can be greatly simpli- 
fied. The differential cross section as a function 
of neutrino energy E^, and struck proton recoil 
kinetic energy Tp (and mass Mp) is 
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In this equation, we have taken {E^ — TpY — E^, 
(i.e., keeping only the lowest order in E^/Mp, 
a very good approximation); the full expres- 
sion was used in the calculations below. The 
neutral-current coupling constants between the 
exchanged Z° and the proton are 

4sin^0,„ 



1 



cv = 



CA = 



= 0.04, 



1.27 



(3) 



where the factor 1.27 is determined by neutron 
beta decay. The cross section for antincutrinos 
is obtained by the substitution ca —ca- As 
will be emphasized below, our results are totally 
independent of oscillations among active flavors, 
as this is a neutral-current reaction. 

We consider only free proton targets; the small 
yield from bound protons creates a small back- 
ground signal, as discussed in Ref. [Q. We use 
the struck proton kinetic energy in the laboratory 
frame as our kinematic variable, as appropriate to 
scintillator detectors. For a neutrino energy E^, 
Tp ranges between and T^"^, where 



The maximum is obtained when the neutrino re 
coils backwards with its original momentum E^ 
and thus the proton goes forward with momen 
tum 2Ey. Since C/i 3> cy, the largest proton re 
coils are favored, which is optimal for detection. 



In this talk, I characterize the supernova neu- 
trino signal in a very simple way, though consis- 
tently with numerical supernova models The 
change in gravitational binding energy between 
the initial stellar core and the final proto-neutron 
star is about 3 x 10^'^ ergs, about 99% of which 
is carried off by all flavors of neutrinos and an- 
tincutrinos over about 10 s. The emission time 
is much longer than the light-crossing time of 
the proto-neutron star because the neutrinos are 
trapped and must diffuse out, eventually escaping 
with approximately Fermi-Dirac spectra charac- 
teristic of the surface of last scattering. In the 
usual model, z/^, Vt and their antiparticles are 
emitted with temperature T ~ 8 MeV, has 
T ~ 5 MeV, and Ve has T ~ 3.5 MeV. The tem- 
peratures differ from each other because Ue and 
have charged-current opacities (in addition to the 
neutral-current opacities common to all flavors), 
and because the proto-neutron star has more neu- 
trons than protons. It is generally assumed that 
each of the six types of neutrino and antineutrino 
carries away about 1/6 of the total binding en- 
ergy, though this has an uncertainty of at least 
50% §. The supernova rate in our Galaxy is es- 
timated to be (3± 1) per century (this is reviewed 
in Ref. 0). 

The expected number of events (assuming a hy- 
drogen to carbon ratio of 2 : 1) is 



N = 70.1 



E 



1 MeV 



T 



10 kpc 
D 



M 



D 



1 kton 



IO- 



CS) 



(Though written slightly differently, this is equiv- 
alent to the similar expression in Ref. ^.) We 
assume D — \Q kpc, and a detector fiducial mass 
of 1 kton for KamLAND. As written, Eq. (|) 
is for the yield per flavor, assuming that each 
carries away a portion E of the total binding 
energy (nominally, Eb = 3 x 10^^ ergs, and 
E — Eb/G)- The thermally-averaged cross sec- 
tion (the integral of the cross section with nor- 
malized Fermi-Dirac distribution) is defined for 
each CH2 "molecule" , and a factor of 2 must be 
included for electron or free proton targets. 



Prior to Ref. [Q, the largest expected yield in 
any oil or water detector was from Pe+P ^ e++n. 
The total cross sections for charged-current + 
p + n and neutral-current ;y + p — > v + p 

have similar forms, though the latter is about 4 
times smaller. However, this is compensated in 
the yield by the contributions of all six flavors, 
as well as the higher temperature assumed for 
and [T = 8 MeV instead of 5 MeV). Thus, the 
total yield from v + p —^ v + p is larger than that 
from i>e+P ^ e++n, when the detector threshold 
is neglected. 

Taking into account radiative, recoil, and weak 
magnetism corrections, the thermally-averaged 
cross section for Pe + P — > e+ -I- n at T = 5 MeV 
is 44 X 10""^^ cm^ (for 2 protons) ||]. These cor- 
rections reduce the thermally-averaged cross sec- 
tion by about 20%, and also correct the relation 
= - 1.3 MeV. The total expected yield 
from this reaction is thus about 310 events in 1 
kton. 

Since the struck protons in v + p ^ v + p have 
a relatively low-energy recoil spectrum, and since 
realistic detectors have thresholds, it is crucial to 
consider the proton spectrum in detail, and not 
just the total yield of neutrinos that interact. 

4. Proton Recoil Spectrum 

The elastically-scattered protons will have ki- 
netic energies of a few MeV. Obviously, these non- 
relativistic protons will be completely invisible in 
any Cerenkov detector like Super-Kamiokande. 
However, such small energy depositions can be 
readily detected in scintillator detectors such as 
KamLAND and Borexino. We first consider the 
true proton spectrum, and then how this spec- 
trum would appear in a realistic detector. 

The proton spectrum (for one flavor of neu- 
trino) is given by 

rl AT fOQ 7 

where f{E^) is a normalized Fermi-Dirac spec- 
trum and the differential cross section is given by 
Eq. (|^). For a given Tp, the minimum required 
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Figure 1. The true proton spectrum in Kam- 
LAND, for a standard supernova at 10 kpc. In 
order of increasing maximum kinetic energy, the 
contributions from i/e, f'e, and the sum of i^^, i^r, 
Vfj^, and Pr are shown with dashed lines. The solid 
line is the sum spectrum for all flavors. Taking 
the detector properties into account substantially 
modifies these results, as shown below. 



neutrino energy is 

[E^jrnin = ^ - V ^ 

The normalization constant C is determined by 
Eq. (||), as the integral of Eq. (^) over all Tp with- 
out the C factor is {a). 

Throughout, I refer to the i^e (T = 3.5 MeV), 
Ve {T — 5 MeV), and the combined t'^, v^, 
and Dr {T = 8 MeV) flavors. Since we know 
that there are neutrino oscillations, this language 
is somewhat incorrect. However, our results are 
totally insensitive to any oscillations among ac- 
tive neutrinos or antineutrinos (since this is a 
neutral-current cross section) , and also to oscilla- 
tions between active neutrinos and antineutrinos 
(since the cross section is dominated by the c\ 
terms). Thus when we refer to the flavor, we 
mean "those neutrinos emitted with a tempera- 
ture T = 3.5 MeV, whatever their flavor compo- 
sition now," etc. The true proton spectra cor- 
responding to the various flavors are shown in 
Fig. |l|. As seen in the flgure, the contributions 



[MeV] 

Figure 2. The quenched energy deposit (equiva- 
lent electron energy) as a function of the proton 
kinetic energy. The KamLAND detector proper- 
ties are assumed. 



of Ve and Ve are quite suppressed relative to the 
sum of v^, Vt, v^, and i>r- 



Figure 3. Analogous to Fig. |l|; the struck pro- 
ton spectrum for the different flavors, but with 
quenching effects taken into account. In order of 
increasing maximum kinetic energy, the contribu- 
tions from i^e, I'ej and the sum of i/^, i^fn and 
Pj- are shown with dashed lines. The solid line is 
the sum spectrum for all flavors. We assume a 1 
kton detector mass for KamLAND. 



5. Quenching 

For highly ionizing particles like low-energy 
protons, the light output is reduced or 
"quenched" relative to the light output for an 
electron depositing the same amount of energy. 
The observable light output Eequiv (i-e., equiv- 
alent to an electron of energy Eequiv) can be 
calculated by integrating Birk's Law with ta- 
bles of dE/dx for protons in the KamLAND oil- 
scintillator mixture The observed energy 

in terms of the proton kinetic energy is shown 
in Fig. 1^. Thus the proton quenching factor 
(Eequzv/Tp) is thus roughly 1/2 at 10 MeV, 1/3 
at 6 MeV, 1/4 at 3 MeV, and so on. 

Using the quenching function shown in Fig. H, 
we can transform the true proton spectrum shown 
in Fig. |l| into the expected measured proton spec- 
trum, shown in Fig. ^. If the quenching fac- 
tor were a constant, it would simply change the 
units of the Tp axis. However, it is nonlinear, 
and reduces the light output of the lowest recoils 
the most. It also reduces the number of events 
above threshold. The anticipated threshold in 



KamLAND is 0.2 MeV electron equivalent en- 
ergy. With the expected proton quenching, this 
corresponds to a threshold on the true proton ki- 
netic energy of 1.2 MeV. The number of events 
above this threshold for each flavor appears in 
Table |^. The measured proton spectrum will pri- 
marily reflect the shape of the underlying Fermi- 
Dirac spectrum for the sum of v^,, i^n, and Dr- 
This has been convolved with both the differen- 
tial cross section (which gives a range of Tp for a 
given E^), and also the effects of quenching. How- 
ever, as we will show, the properties of the initial 
neutrino spectrum can still be reliably deduced. 

Background considerations, while important, 
are a small correction, and so we ignore them 
here; see Ref. |^ for a complete discussion. 

6. Proton Spectrum Fits 

The measured proton spectrum can be used 
to separately determine the total energy of the 
h'^, Vr, f'/i, and i>r neutrinos and their time- 
averaged temperature. The total number of de- 
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Table 1 

Numbers of events in KamLAND (1 kton mass 
assumed) above the noted thresholds for a stan- 
dard supernova at 10 kpc, for the separate flavors 
or their equivalents after oscillations. Oscillations 
do not change the number of neutrinos at a given 
energy, and the neutral-current yields are insensi- 
tive to the neutrino flavor. Equipartition among 
the six flavors is assumed (see the text for discus- 
sion). The thresholds are in electron equivalent 
energy, and correspond to minimum true proton 
kinetic energies of and 1.2 MeV. 
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tected events is proportional to the portion of the 
total binding energy carried away by these four 
flavors, and we denote this by £'*°* (note that 
this is not the total binding energy Eb)- For a 
standard supernova, = 4(£:b/6) = 2/3Eb 
2 X 10^'^ ergs. We denote the temperature of these 
four flavors by T. If only the total yield were mea- 
sured, as for most neutral-current reactions, there 
would be an unresolved degeneracy between 
and T, since 



N ^ E*' 



T 



El\ then (a) 



(8) 



T". For 



Note that for a 
ly + d iz+p+n in SNO, for example, a ~ E'^, so 
N ~ E*°*T. Thus for a given measured number 
of events, one would only be able to define a hy- 
perbola in the plane of and T. The scaling is 
less simple here because of threshold effects, but 
the idea is the same. 

Here we have crucial information on the shape 
of the neutrino spectrum, revealed through the 
proton spectrum. To remind the reader, in most 
neutral-current reactions there is no information 
on the neutrino energy, e.g., one only counts the 
numbers of thermalized neutron captures, or mea- 
sures nuclear gamma rays (the energies of which 



depend only on nuclear level splittings). 

We performed quantitative tests of how well the 
parameters E*°* and T can be determined from 
the measured proton spectrum. (We did also in- 
vestigate the effects of a chemical potential in the 
Fermi-Dirac distribution, but found that it had 
little effect. This is simply because the cross sec- 
tion is not rising quickly enough to see the tail of 
the thermal distribution in detail .) Of course, 
if the distance to the supernova is not known, 
then we are effectively fitting for E*°* /D^ . 

We performed Monte Carlo simulations of the 
supernova signal in KamLAND and made chi- 
squared fits to determine and T for each fake 
supernova. To perform the fits, we started with 
an "ideal" spectrum, as described by the integral: 



fd^r 



c 



dT;G{T;;Tp,STp) 



da 



X / dE,f{E,)—{E,,T;) (9) 



where the inner integral is as in Eq. (^, with 
the addition that quenching corrections are ap- 
phed to Tp after convolution with f{E^). For 
the Gaussian energy resolution GiTp, Tp, STp), we 
used STp = 0.1y/Tp/{l MeV) 0. The nor- 
malization constant C is given by comparison to 
Eq. (H). Example spectra are shown in Fig. |[ 

Using {dN/dTp)ificai, we binned the spectrum 
by the following integral: 



N.= r"^^^^^dTp(^\ (10) 

JiE^,^), ideal 

where Ni is the number of events in bin i, and 
{Emin)i and {Emax)i are the minimum and max- 
imum energies for bin i. Eight bins of variable 
width were used, chosen to contain roughly the 
same number of expected events per bin. For a 
chosen E*°* and T, this was the starting point of 
our Monte Carlo (and the bin boundaries were 
kept fixed) . For each fake supernova, we sampled 
the number of events in each of these bins ac- 
cording to the appropriate Poisson distributions. 
The resulting spectrum was as one might obtain 
from a single supernova, given the finite number 
of events expected. We then varied and T 
in Eq. (nn) until the values that best fit the fake 
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Figure 4. Example spectra with different values 
of and T, all chosen to give the same num- 
ber of events above an electron equivalent thresh- 
old of 0.2 MeV (true proton energy 1.2 MeV) in 
KamLAND. Though not shown in this figure, the 
spectrum above 2 MeV is included in our analy- 
sis. At the 0.2 MeV point, from left to right these 
correspond to {E^°\T) = (4.2, 6), (2.0, 8), (1.4, 
10), respectively, with in 10^^ ergs and T in 
MeV. 



spectral data were determined. For a given set 
of assumed and T, this procedure was re- 
peated many times. The distributions of the final 
and T thus reveal the expected errors on fit- 
ting E*°* and T for a single real future supernova. 
Three examples are shown in Fig. |^, where one 
can see that and T can each be determined 
with roughly 10% error. 

7. Discussion and Conclusions 

We have shown that neutrino-proton elastic 
scattering, previously unrecognized as a useful 
detection reaction for low-energy neutrinos, in 
fact has a yield for a supernova comparable to 
Ve + p ^ + even after taking into account 
the quenching of the proton scintillation light and 
assuming a realistic detector threshold. 

In addition, the measured proton spectrum is 
related to the incident neutrino spectrum. We 
have shown explicitly that one can separately 
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Figure 5. Scatterplot of 10"^ fitted values, in the 
and T plane, for the labeled "true" values, 
where £'*°* is the total portion of the binding en- 
ergy carried away by the sum of v^^ D^, and 
Dr, and T is their temperature. The values of 
and T were chosen such that the numbers of 
events above threshold were the same. The mea- 
sured shape of the proton spectrum breaks the de- 
generacy between these two parameters. Without 
that spectral information, one could not distin- 
guish between combinations of and T along 
the band in this plane that our three example re- 
gions lie along. 



measure the total energy and temperature of f^, 
Vr, Vfj,, and Pt, each with uncertainty of order 10% 
in KamLAND. This greatly enhances the impor- 
tance of detectors like KamLAND and Borexino 
for detecting supernova neutrinos. 

For Borexino, the useful volume for supernova 
neutrinos is 0.3 kton, and the hydrogen to car- 
bon ratio in the pure pseudocumene (C9H12) is 
1.3 



1 |12 



so there are about 4.7 times fewer 
free proton targets than assumed for KamLAND. 
However, the quenching is less in pure scintilla- 
tor (KamLAND is about 20% pseudocumene and 
80% paraffin oil ) , and the errors on £'*°* and 
T scale as 1/ \/N , so that the precision in Borex- 
ino should be about 20% or better. 

Other techniques for bolometric measurements 
of supernova neutrino fluxes have been stud- 
ied. Detectors for elastic neutral-current neu- 
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trino scattering on electrons [|l3| and coherently 
on whole nuclei have been discussed, but 
never built. If neutrino oscillations are effec- 
tive in swapping spectra, then the temperature 
of the "hot" flavors may be revealed in the mea- 
sured positron spectrum from Ve + p ^ + n; 
two recent studies have shown very good preci- 
sion (< 5%) for measuring the temperatures and 
the total binding energy However, they 

assumed exact energy equipartition among the 
six neutrino flavors, whereas the uncertainty on 
equipartition is at least 50% Nevertheless, 
under less restrictive assumptions, this technique 
may play a complementary role. Finally, since for 
different cross sections, the neutral-current yields 
depend differently on temperature, comparison 
of the yields may provide some information . 
However, there are caveats. In neutrino-electron 
scattering, the neutrino energy is not measured 
because the neutrino-electron angle is much less 
than the angular resolution due to multiple scat- 
tering. The scattered electrons, even those in a 
forward cone, sit on a much larger background of 
i>g + p + n events, so it is difficult to mea- 

sure their spectrum also, their total yield 

is only weakly dependent on temperature. At 
the other extreme (see Fig. 3 of Ref. jl^), the 
yield of neutral-current events on ^^O de- 
pends strongly on a possible chemical potential 
term in the thermal distribution. 

It is important to note that the detection of 
recoil protons from neutron-'pTotoii elastic scat- 
tering at several MeV has been routinely ac- 
complished in scintillator detectors (see, e.g., 
Ref. iQ). Since both particles are massive, the 
proton will typically take half of the neutron en- 
ergy. This reaction provides protons in the same 
energy range as those struck in neutrino-proton 
elastic scattering with E,y ~ 30 MeV. This is a 
very important proof of concept for all aspects of 
the detection of low-energy protons. 

Though low-energy backgrounds will be chal- 
lenging, it is also important to note that the back- 
ground requirements for detecting the supernova 
signal are approximately 3 orders of magnitude 
less stringent than those required for detecting 
solar neutrinos in the same energy range (taking 
quenching into account for our signal). Borcxino 



has been designed to detect very low-energy so- 
lar neutrinos, and KamLAND hopes to do so in 
a later phase of the experiment. 

These measurements would be considered in 
combination with similar measurements for 
and De from charged-current reactions in other 
detectors. Separate measurements of the total 
energy and temperature for each flavor will be 
invaluable for comparing to numerical supernova 
models (|J2^ . They will also be required to make 
model-indcpcndcnt studies of the effects of neu- 
trino oscillations |2|] . If the total energy release 
Eb in all flavors has been measured, then 
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thus allowing a direct and unique measurement of 
the newly-formed neutron star properties, princi- 



pally the mass M^vs |23 
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